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Tall buildings with irregular shapes and considerable heights are gaining popularity 
in creating the vertical cities around the world. They also considered one of the 
major energy consumers with little regards to sustainability. Tall building is wind- 
sensitive structure and shape plays major role in determining wind loads, which 
usually govern the design of its lateral resisting system. Thus, evaluating wind loads 
properly and designing an optimal lateral system accordingly are the main 
challenges attributed to the design process of irregular tall building. This paper 
presents a computational procedure for the optimal design of wind-resistant 
irregular tall building to minimize the total weight of structure within design 
requirements in single digital environment. That is achieved firstly by creating a 
digital system of computational fluid dynamic (CFD) analysis that is coupled with 
pressure-load translation (PLT) algorithm to evaluate the wind motions on irregular 
tall buildings and generate the design wind loads accordingly. Genetic Algorithm 
(GA) with enhance design constrains function of lateral displacements, inter-story 
drifts and top acceleration is then developed to perform structural optimization. A 
numerical example using 70-story twisting reinforced concrete building is 
implemented to verify the feasibility of the developed computational 
procedures. Steady and incompressible flow applied at (0°) angle of attack was 
implemented in the CFD model to simulate the wind flow on the studied building. 
Genetic algorithm with improved design constraints of static and dynamic design 
requirements was developed to optimize the structure effectively and efficiently. 
The numerical example demonstrates its effectiveness by achieving 35.71% 
reduction of concrete volume from the original lateral structural system design. 
This is also translated into a sustainability value by lessening the embedded carbon 
dioxide by 4,400 tons. 
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1 Introduction 


Social and economic accelerating development, together with new 
engineering and construction techniques have inspired architects and 
engineers to design atypical, complex and slender tall buildings to 
reflect their creativity and innovation (Alkhatib et al., 2021a). These 
developments are associated with complexity in wind-induced 
motions and its attributed structural factors such as low natural 
frequency, lack of sufficient damping and increment in flexibility and 
slenderness (Sharma et al., 2018), where wind-induced forces usually 
govern the design of tall building (Elshaer et al., 2017). Thus, finding a 
balance between safety, serviceability and economy has been a 
challenge faced by architects and engineers when designing 
modern irregular tall buildings. This challenge is leveled up when 
sustainability factor is added into the equation by producing energy- 
efficient design. 

Because building shapes plays critical role in deriving its 
aerodynamic characters, modification of tall building forms 
proved to be effective in mitigating wind-induced loads on tall 
buildings (Xie, 2014; Elshaer et al., 2017). Several studies 
investigated different aerodynamic modifications on tall 
buildings including major modifications such as tapering, 
twisting and setback (Kim and You, 2002; Kim et al., 2008; 
Xie, 2014; Bairagi and Dalui, 2018) that have a significant impact 
architecturally and structurally on building design, or minor 
modifications such as corner configuration (Zhengwei et al., 
2012; Li et al., 2018; Thordal et al., 2020; Mandal et al., 2021), 
where it has little to no effect on the intent design, yet still shows 
an effective results in reducing wind responses. While Kim and 
You (Kim and You, 2002) and Xie (Xie, 2014) found that tapering 
is more effective in mitigating across-wind responses, Xie (Xie, 
2014) further 
demonstrating that it has and extra advantage of equalizing 


explored twisting shape modification 
winds distribution and responses along wind directions. 
Minor corner modifications, including rounding, chamfering 
and recession, were also investigated by Tamura et al. 
(Tamura et al., 1998) and Elshaer et al. (Elshaer et al., 2017) 
and founded to be effective in reducing wind responses. Li et al. 
(Li et al., 2018) demonstrated that chamfered corners resulted in 
an effective reduction of a long-wind loads and recessed corners 
showed similar effectiveness in reducing across-wind responses, 
but Mandal et al. (Mandal et al., 2021) proved that rounded 
corners are more effective in reducing wind loads than chamfered 
ones. Alkhatib et al. (Alkhatib et al., 2022a) introduced a new 
computational method to optimize wind loads on irregular tall 
buildings with minor aerodynamic modifications to the external 
shape of the buildings. Touloupaki et al. (Atabay, 2009), in the 
other hands, proposed a parametrical approach to minimize 


energy consumption of tall building by altering its geometrical 
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form. Hassanli et al. (Hassanli et al., 2018) and Ruiz et al. (Ruiz 
et al., 2021) investigated the effectiveness of tall building 
aerodynamic treatements to wind energy harvesting by 
incorporating wind power system directly into the built 
environment. Although these studies present a valuable 
insight by exploiting computational methods to evaluate wind 
tall buildings properly, 
configuration is not architecturally and structurally preferred. 


motions on altering building 

Quite the opposite, structural optimization does not rely on the 
changing of the building external geometries or requires additional 
cost. However, most of previous research work of structural 
optimization focused on seismic loading, with little regard to 
wind loading (Huang et al., 2015; Fu et al., 2018; Alkhatib et al., 
2022b). In that context, Chan et al. (Chan and Chui, 2006; Chan 
et al., 2009) developed a set of structural optimization methods based 
on Optimally Criteria (OC) for optimal design of tall buildings 
subjected to wind actions and constraints. Huang (Ming-feng, 2013) 
as well, proposed a performance-based design approach by 
employing OC for developing optimal structural design of tall 
building satisfying wind-resistance requirements. However, Fu 
(Fu et al., 2018) demonstrated that employing OC approach may 
lead to inaccurate result of frequency constraint because strain 
energy does not account for the gradient of internal forces to its 
design variables and proposed an eigen value approach for 
computing the sensitivity of frequency constraints. Moreover, OC 
requires wind-induced constrains to be explicitly expressed in 
complicated mathematical formulas, making it not practically 
applicable for industrial application. Additionally, dealing with 
complex and large-scale models, structural response is naturally a 
non-linear function to the design variables and OC may not lead to 
discrete solutions and may trapped into local optima. It is also 
observed that most of the optimization techniques used today for 
solving engineering problems accept the design variables to be 
“continuous” and such characteristics cannot be employed in all 
engineering problems. For architectural problems, besides the 
mathematical formulation techniques, using “discrete” variables 
based algorithm is often deemed necessary (Atabay, 2009; Abd 
Elrehim et al., 2019). 

Taking these interesting studies as a benchmark for this 
research to address their shortcoming, this study presents a 
new integrated approach for wind-induced structural 
optimization of irregular tall building within single digital 
environment. In the first place, a computational fluid dynamic 
(CFD) system coupled with pressure-load translation (PLT) 
algorithm is developed to evaluate wind motions on irregular 
tall building and generate the corresponding design wind 
loads. Next in order, structural optimization strategy is 
developed using Genetic Algorithm, with its ability for 
global search and efficient parallelism, to find the optimal 
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weight of the structure. More importantly, because lateral 
drifts and top acceleration are the focus in the design of wind- 
resistance tall building, an enhanced design constrains penalty 
function is developed and integrated within GA process. Static 
and dynamic characteristics of the building are then assessed 
through finite element (FE) integrated solver to evaluate 
structural responses against design constraints pre-defined 
limits. 

The paper is structured into 5 sections. Section 1 (this 
section) describes the topic under investigation and 
summarizes the related studies with its shortcomings to be 
addressed in this work. Section 2 demonstrates CFD wind 
loads generation and improved GA structural optimization 
methodologies. Section 3 verifies the proposed approaches by 
implementing a numerical example on real existing irregular 
tall building where analysis and results are presented and 
discussed. Finally, Section 5 draws the conclusion of this study 
and highlights the future recommendation for extending this 
research work. 


2 Optimization methodology 
2.1 Computational wind loads evaluation 


Computational Fluid Dynamic (CFD) has recently 
extended its applications to architectural problems 
including tall building aerodynamic studies (Mou et al., 
2017; Paul and Dalui, 2021; Sanyal and Dalui, 2022) and 
this opens an opportunity to utilize CFD in the preliminary 
stage to evaluate wind motions on irregular tall buildings. 
However, the idea in this work is to develop an integrated 
optimization approach where CFD output can be translated 
into static wind loads acting on each story of a given structure 
within a single environment in an automated manner. Thus, 
parametric modelling using Rhino3D (Rutten and McNeel, 
2007) and its embedded generative design tool Grasshopper 
(GH) is employed firstly to generate the initial model 
representing the bluff body of building exterior shape. 
OpenFOAM (Weller et al., 1998), an open-source 
computational fluid dynamic solver that is widely and 
reliably used in academia and industry to solve 
computational problems, is utilized to perform the outdoor 
airflow analysis. However, to run OpenFOAM CED cases 
within GH environment, a python wrapper for OpenFoam 
C++ library is used. Significant parameters to set up CFD 
environment are defined and evaluated to ensure proper 
evaluation of wind pressure on the selected building. That 
includes turbulence models, mesh resolutions, wind profile, 
wind speed and convergence criteria. 

Pressure-Load Translation (PLT) algorithm is then 
developed to translate CFD output pressures into static 
wind loads. The PLT technique starts by transforming the 
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pressure exerted by the building exterior façade mesh surfaces 
into vector forces with respect to its normal directions. Sorting 
algorithm is written by Python (vanRossum, 1995) to combine 
and group the collected vectors horizontally and vertically for 
each story height that is measured by mesh divisions. For each 
story, two vectors are generated to represent along-wind force 
(Fy) and across-wind forces (Fx) with respect to building 
major axis directions, while torsional moments (Mz) are 
computed as the output of their cross-product to the 
building’s centroid. Generated Forces and moments are 
then assigned to the diaphragm center of the structure for 
each respective story. The overall process of CFD optimal 
wind loads generation is shown in Figure 1. 


2.2 GA structural optimization 
methodology 


Because satisfying serviceability criteria become more 
challenging than achieving strength requirement (Chan 
et al., 2009), also since overall design of tall building is 
usually governed by lateral stiffness, the output sizes from 
lateral stiffness serviceability-based optimization are larger 
than strength based design sizes (Soegiarso and Adeli, 1997; 
Pezeshk et al., 2000). Hence, only serviceability requirements 
are employed as constraints for this optimization study, where 
strength constraints can be dealt with independently (Chan 
et al., 1995). Practically, strength and capacity design are 
regularly advised to be performed on a member-by-member 
approach, with reference to the standard code adopted (Zou 
and Chan, 2005). 


2.2.1 Objective function 

In structural design optimization, generally weight of 
structure is assumed to be the objective function. 
Considering concrete tall building structure, the objective 
function can be expressed as minimization of the structural 
lateral system concrete weigh. Thus, for a 3D structural frame 
consisting of I = 1,2,3 ..., N types of core walls, columns and 
beams, and j = 1,2,3....M stories, the wind-resistance 
objective function is formulated as: 


To minimize: 
N 
f(x) = YpAil; 
şil 


Limited to: 


b! <bi<bY; ht ch sh? st! stit? 


i = t= 


bi, hi, ti = S; € {S} 


Where p is concrete density, A; is sectional area for the ith 
type of element section, and J; is the length for ith type of element 
section. b;, h; are the width and depth for column/beam section in 
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FIGURE 1 
Computational wind loads generation process 


ith type section respectively and t; is the thickness of corewall in 
ith type of section. S; is the indecies for the section sizes domain 
defined within {S}. 


2.2.1.1 Design constraints 

General formulation of constrains function for k = 1,2,3 
...,V number of constrains, for j = 1,2,3....,M number of 
stories is expressed as follow: 


9x (x) <0 


Giving the problem of wind-resistance optimization, 
lateral drifts and top acceleration are expressed as the main 


design constraints restricting their values within 
recommended limits. 
UH vH Vui + vi 
= < Ĝlimit 5 < limit ; < Olimit 
H H H 
Uj — Uj- 
j 7 Uja 
< Aiimit , 
Zj = Zj-1 
y 2 
yy uj- uji) +(vj- vj) 
j j-l i ( J J J d 
< Alimit , < Aiimit 
Zj Zj-1 Zj = Zj-1 
XH Š Alimit 3 Yy S Amit} XV y S limit 
However, GA can only work with unconstrained 


optimization problem, thus, constrains were transformed into 
normalized form: 
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S| 9 
= Z. 
a) 6 
S =] = 
sl] S| & 
ST | =) 
fl] Sl x 
un vH mar 
a -1<0; *--1<0; 5 -1<0 
limit limit limit 
ujuja YjYja (ujuja) + (vzv) 
Zj-Zj-1 Zj-Zj-1 Zį-Zj-1 
——-1<0; == -1<0; L -1<0 
Altimit Aiimit Aiimit 
om = z 
-1<0; YH _1 <9, H_1 <0 
limit limit limit 


uy, vH are the along-wind and across-wind maximum 
horizontal displacement at top of the structure and their 
diagonal vector, djimit is the allowable limit stipulated by 
the code (Eurocode H/500), uj, vjare the along-wind and 
across wind inter-story drift and their diagonal vector at 
jth floor in respect to floor height zj, Ajimiy is the allowable 
inter-story drift. xy, Yy and Armi are the top story peak 
acceleration in x and y directions and their recommended 
limit in a adopted standard, respectively. Threshold of top 
peak acceleration is computed by employing the lower natural 
frequency generated by eigenvalue analysis as per ISO 6897 (I 
S O., 1984) and ISO 10137:2007 (ISO Organization for 
Standardization, 2007): 


Xq = \2In fı T + (08+ =R) exp ( -0.35 — 0.41 1n fı) 


Xu is the peak acceleration for T minutes within a return 
years period R, f is the natural frequency of first mode, 
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FIGURE 2 
Genetic Algorithm optimization process. 
2In fı T is peak factor, (0.68 +28) is adjustment factor 10 years. XVq is adopted to evaluate the complex motion of 
and exp (-0.35-0.41lnf;) is natural frequency irregular tall building to estimate its corresponding horizontal 
acceleration function root mean square. Usually, acceleration as recommended by (Melbourne and Palmer, 
observation time T is 10 min and return period R is 1, 5 or 1992). 
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2.2.1.2 Fitness function 
In order to calculate fitness value for each individual, and 
because constrains were normalized, the objective function has to 
be modified by a violence coefficient, C: 
g(x) if gx(x)>0 
Ck = 
0 if ge(x)<0 


Where C = pcx 
The modified objective function can then be written as: 
@(x) = f(x) (1+ KC) 


Where ®(x) is the modified objective function, f (x) is the 
actual objective function, C is the violation coefficient and K is 
a penalty factor with a value depending on the problem. In 
this case study, K is assumed to be 10 for deflection constrain 
and 1,000 for drift constrain. The objective function is 
then required to be transformed into a fitness value. 
(Goldberg, 1989) 
minimization problem, to generate fitness value, ® (x) shall 


According to for optimization 
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be subtracted from large fixed value, so that all fitness values 
will be positive, and individuals with minimum criteria will 
have higher fitness values. The fitness value can then be 
obtained by subtracting the value of unconstrained 
objective function ®; (x) from the summation of maximum 
and minimum unconstrained objective function values within 
a population consisting of all individuals in ith generation as 


follow: 


F; = [®; (x) max + QD; (x) min ] g Q; (x) 


Fitness factor then can be derived as: 


Fe; = F; f Fag. 
Fay. =) Fifn 


where n is total number of individuals in a population. The 
overall process of genetic algorithm structural optimization 
procedures is illustrated in Figure 2. 


floor plans 3 -21 


FIGURE 3 
Cayan tower dubai (Alkhatib et al., 2022a). 
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FIGURE 4 
Parametric models: (A) Bluff, (B) detailed and (C) structural. 


FIGURE 5 
Computational domain (outer) meshing. 


3 Numerical example 


Case Study: Cayan tower in Dubai, as shown in Figure 3, the 
tallest twisting tower in the world, present a good example of the 
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FIGURE 6 
Along-wind forces (Fy) along y-direction. 


current trend of complexity in tall building design. The tower is 306 m 
in height, with spiral twisting shape of 90° along its height, inspired by 
human DNA and designed by SOM architectural group (Golasz- 
Szolomicka and Szolomicki, 2019). The lateral structure system of this 
concrete building consists of central core wall and perimeter moment 
frame. To simplify the studied problem, typical story height of 4 m is 
assumed for all 75 story with total structure height of 300 m. 


3.1 Wind loads generation 


Two versions within the same parametric model are created as 
illustrated in Figure 4 (a) and (b). The first one is a bluff model 
representing the exterior shape of the building to be used for CFD 
wind analysis stage. The other one is element detailed model to be 
translated into structural model, Figure 4C to be integrated into FE 
solver for structural responses evaluation during optimization stage. 

To set up the CFD environment, different parameters and 
settings adopted to ensure proper analysis and reliable output. 
In the first place, domain size of virtual tunnel, shown in 
Figure 5, was determined based on building height (H) with 
2.3 H, 10 H, 2.3 H and 2.3 H for windward, leeward, sides and 
top, respectively as recommended by (Franke et al., 2007). 
Blockage ratio was then computed at 1.06%, below the 
threshold ratio of 5%, to provide sufficient space for the air 
flow to be fully developed (Mou et al., 2017). The Outflow 
boundary condition was employed to simulate steady and 
incompressible wind flow. Wall boundary conditions were 
used to describe the wall of the computational domain and the 
surface of the studied building. SIMPLE pressure-velocity 


frontiersin.org 


Alkhatib et al. 


10.3389/fenrg.2022.1017813 


Building Story 


Fx (KN) 


FIGURE 7 
Across-wind forces (Fx) along x-direction. 


Building Story 


0 100 200 300 400 500 600 
Mz (kN.m) 


FIGURE 8 
Torsional moment (Mz) around the building's center. 


coupling algorithms was adopted for solver setting, and 
Second-Order Upwind Scheme was adopted to discretize 
momentum, turbulent kinetic energy, and turbulence 
dissipation rate (Meng et al., 2018). Finally, convergence of 
simulation is assumed to be reached when residuals of x, y, z 
momentum, k, e and continuity are less than 10° °°. 
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FIGURE 9 
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FIGURE 10 
Crossover mechanism. 


FIGURE 11 
Mutation mechanism. 


Realizable k—e turbulence model was adopted to perform 
numerical simulation. That is because it is sensitive to rapid 
strain and streamline curvature, flow separation, re-attachment 
and recirculation (Ozmen et al., 2016), while at the same time, it 
allows mathematical constrains to fit Reynolds stresses (Tutar 
and Oguz, 2002). For mesh resolutions, maximum grid size is 
determined by dividing shortest width of the building over 
10 and multiplying the output result by V2 for each refinement 
needed (Alkhatib et al., 2022a). Two staged refinements were 
done based on similar study (Alkhatib et al., 2022a) as no 
considerable change of output values is observed for advanced 
refinement, while reasonable computational time is maintained. 
The selected mesh size is computed by 1.54 m reaching to a total 
of 972,660 cell numbers within the computational domain. 
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TABLE 1 Structural sections grouping. 
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Structural element Group Original size Lower limit Upper limit 
Core Wall Core wall thickness (mm) 
CW1 350 725 1,500 
Cw2 »250 625 1,400 
CW3 200 575 1,350 
CW4 »150 525 1,300 
CW5 100 475 1,250 
Column Column Section (Depth x Width) mm 
CL1 500 x 1,000 875 x 600 1,650 x 1,100 
CL2 475 x 1,000 850 x 575 1,625 x 1,100 
CL3 450 x 975 825 x 550 1,600 x 1,075 
CL4 425 x 950 800 x 550 1,575 x 1,050 
CLS 400 x 950 775 x 525 3950 x 1,050 
Beam Beam Section (Depth x Width) mm 
BM1 500 x 1,000 875 x 600 3650 x 1,100 
BM2 450 x 975 825 x 550 :600 x 1,075 
BM3 400 x 950 775 x 525 3950 x 1,050 
BM4 350 x 900 725 x 500 3900 x 1,000 
BM5 »300 x 875 675 x 450 :450 x 975 


The generated wind loads for the typical 0° wind direction 
that was adopted in this research work are shown in Figure 6, 
Figure 7 and Figure 8 for along-wind forces (Fy), across-wind 
forces (Fx) and torsional moment (Mz), respectively. It can be 
observed that greater forces are generated along the angle of 
attack in comparison to the across-wind forces. Torsional 
moment experienced some disruptions due to the nature of 
the building with its irregular verticality especially at the mid- 
height of the building. 


3.2 GA structural optimization 


To implement the developed genetic algorithm technique in 
the studied example, several configuration and parameters 
should be established. That includes problem formulation of 
the investigated problem to suite the optimization algorithm, 
coding and decoding strategy and controlling parameters. 


3.2.1 Problem formulation 

The investigated structure was firstly grouped and arranged 
to simplify the optimization process. Structure was divided into 
5 groups along its height. Group 1, group 2, group 3, group 
4 and group 5 represent structural elements from 1st to 15th 
floor, 16th to 30th floor, 31st to 45th floor, 46th to 60th floor 
and 6lst to 75th floor, respectively as illustrated in Table 1. 
Since lateral structural system is the aim for this optimization, 
other structural elements including slab and internal columns 
are assumed to be constant throughout the height of the 
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structure. Original structural elements are set based on 
reference structural plan of the case study from (Golasz- 
Szolomicka and Szolomicki, 2019) and (Marmar et al., 2021) 
together with engineering experience and judgement. Variation 
for lower and upper limits are set to be within Eurocode design 
guidance and standard (Hendy and Smith, 2007). To alleviate 
long computational processing time by limiting number of 
design variables, columns and beams are assumed with fixed 
ratio of 3:2 for depth and width respectively. Thus, sections’ 
width is changed in proportion to their depth by the fixed ratio. 
Finally, taking constructability into account, the variance in 
sections sizes is determined by 25 mm intervals. 

Design constraints limitations for structural displacement 
and inter-story drift are adopted based on Eurocode Standard 
(Hendy and Smith, 2007). The computed limits are 600 mm and 
0.0015 for structural displacement and inter-story drifts, 
respectively. Top acceleration limit of 0.21 m’/s is generated 
based ISO 10137:2007 
Standardization, 2007) with 1 year return period. 


on (ISO Organization for 


3.2.2 Coding and decoding 

Several ways to represent design variables in GA including 
real and binary values. This research however adopted binary 
method, as shown in Figure 9, for its highly efficiency uniform 
crossover (Hu and Di Paolo, 2008), and faster computing time as 
it deals with machine original language of 0 and 1. For the 
15 design variables represented by the assigned groups, change of 
design variables sectional area were planned in a set of pool with 
maximum of 625mm, 150mm in section size reduction or 
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increment respectively, and within 25 mm intervals. This will 
produce a set {S} of 32 sections variance: 


{S} = {625, 600, 575, . . . , —125, —150} 


To produce a binary code for {S}, 5 bits is needed for 
32 binary coding with (00000) representing first element in 
set {S}and (11111) representing the last element in set {S}. For 
the 15 design variables (groups), each of them is represented by a 
gene (substring) corresponding to a value from the {S} set and its 
binary coding. All multi design variables genes are then linked to 
form chromosome (string). Therefore, for the given 15 design 
variables with each substring consisting of 5-bits, each individual 
solution represented by a chromosome (string) having a length of 
75-bits of 0’s and 1’s. 


3.2.3 GA control parameters 

The main controlling parameters for genetic algorithm are 
population size and generations number. Population size 
projects number of individuals in each generation, and 
number of generations projects maximum number of 
iterations for computational process of genetic algorithm, 
unless convergence is achieved or stopping criteria is 
activated. Considering the problem given in this case study, 
size of design variables, computational time and pre-test of the 
of 
generations were set to be 50 each. After initial population is 


developed algorithm, population size and number 
randomly created and evaluated, the next generation is 
produced by undergoing through elitism, crossover and 
mutation process, as shown in Figures 10, 11. 

Elitism was adopted in the optimization process because of 
its effectiveness in improving the efficiency of GA (De Jong, 
1975). Elitism is achieved by preserving the highest performance 
individual solution from one generation to another by simply 
copying the fittest chromosome to the next generation before all 
individuals in the population undergo crossover and mutation 
process. 

Crossover methods including single-point, multi-point, 
uniform and variable-to-variable crossover were tested to find 
the most effective method for the given optimization problem 
(Hasangebi and Erbatur, 2000). Multi-point crossover was 
adopted giving that it has a better logic for the optimization 
problem in hand. In multi-point crossover, parents’ 
of 
substrings at random locations based on the introduced 
(cp). 
interchanged between parents to produce new offspring. That 


chromosomes (strings) are decomposed to number 


cutting points Substrings are then swapped and 
allows the algorithm to produce new combination of structural 
sections along structure groups to be evaluated. Crossover rate 
which decides the possibility of crossover to happen in a given 
generation is set to be 0.7 as generally implemented in GA 
optimization. 

Mutation was also applied to introduce and maintain 
diversity in exploration of new solutions. Bit-flip is the 
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method applied in binary coding problems and it is 
achieved by randomly selecting one individual chromosome 
from in a given population and randomly flipping one of its 
gene’s values (from 0 to 1 and vice versa). Mutation rate which 
decides the possibility of mutation to happen in a given 
generation is set to be 0.17 as generally implemented in GA 
optimization. 


3.3 Analysis and discussion 


Design optimization cycle started by generating the random 
population and undergoing evaluation process by implementing 
the above-mentioned objective function and its design contains 
That is 
serviceability analysis and dynamic eigen value analysis. The 


mathematical formulations. involved both static 
first critical modes of shape are shown in Figure 12. It is observed 
that due to the complexity of the building form, critical mode of 
shapes also shows similar complexity in deformation, which 
cannot be predicted normally as regular shapes. This proves 
the necessity of the developed formulation to take into account 
the deformation in both x and y direction, as well as in their 
diagonal resulted vector direction. 

For each generation, the fittest individual within the 
population is ranked first and documented. That is to allow 
its use for the second generation if there is no better solution is 


FIGURE 12 
Eigen-value modal analysis (critical mode of shapes). 
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TABLE 2 Genetic Algorithm optimization summary. 


Generation Selected Constraints along X-direction Constraints along Y-direction Constraints along XY-direction Objective Penalty Fitness 
number individual value index 
Deflection Drift Acceleration Deflection Drift Acceleration Deflection Drift Acceleration 
(mm) (m2/s) (mm) (m2/s) (mm) (m2/s) 
Original Original 242.25 0.00147 0.119 387.00 0.00182 0.115 456.55 0.00234 0.166 30984.13 0 Original 
1 39 279.13 0.00174 0.119 473.47 0.00229 0.115 549.63 0.00287 0.166 23775.69 0 1.020 
2 37 292.87 0.00184 0.119 490.04 0.00236 0.115 570.89 0.00299 0.166 23485.89 0 1.154 
3 22 275.85 0.00173 0.119 479.72 0.00227 0.116 553.37 0.00285 0.166 23233.10 0 1.152 
5 7 277.02 0.00171 0.120 490.26 0.00227 0.117 563.11 0.00284 0.167 22580.85 0 1.170 
6 8 304.74 0.00190 0.121 513.58 0.00242 0.117 597.18 0.00307 0.168 21338.79 0 1.020 
15 45 303.35 0.00187 0.121 512.59 0.00239 0.117 595.63 0.00303 0.168 21052.61 0 1.020 
22 19 291.13 0.00180 0.120 502.33 0.00233 0.117 580.60 0.00294 0.168 20724.66 0 1.042 
24 31 299.22 0.00184 0.121 512.73 0.00237 0.117 593.65 0.00300 0.168 20602.64 0 1.020 
28 46 297.01 0.00185 0.120 508.69 0.00238 0.116 589.05 0.00302 0.167 20541.63 0 1.020 
33 30 300.75 0.00188 0.120 512.98 0.00240 0.116 594.64 0.00305 0.167 20465.36 0 1.165 
39 43 299.76 0.00187 0.120 511.11 0.00240 0.116 592.53 0.00304 0.167 20395.15 0 1.158 
41 14 302.33 0.00189 0.120 514.54 0.00241 0.117 596.79 0.00306 0.167 20292.68 0 1.062 
43 9 302.16 0.00189 0.120 515.12 0.00241 0.116 597.20 0.00306 0.167 20192.11 0 1.063 
50 10 302.58 0.00189 0.120 516.36 0.00241 0.117 598.48 0.00306 0.167 19917.11 0 1.157 
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FIGURE 13 
Core-wall section optimization. 
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FIGURE 15 
Beam section optimization. 
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FIGURE 14 
Column section optimization. 
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FIGURE 16 
Structural concrete volume optimization. 


Objective Original (m Optimal (m Reduction (m°) Reduction (%) CO, (tons) 
(Elshaer et al., (Elshaer et al., 
2017)) 2017)) 

Lateral Structural System 30,984 19,917 11,076 35.71 4,400 


found, hence algorithm does not lose the best in hand. 
Summary of the optimization cycle and its associated 
objective and constraints values are recorded in Table 2. To 
minimize the lengthy list of the produced 50 generations, and 
for better result representation, only generations that 
experienced a change in their individual solutions output 
are depicted in Table 2. 

Several valuable insights can be discussed from the recorded 
solutions stipulated in Table 2 above. For the design constraints, 
it was observed that structural displacement governs the design 
and optimization, while drift and top acceleration still stand far 
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from their limited values. This provides engineers with valuable 
information on the critical structural response to be dealt with. 
Besides, it also opens an opportunity for extended optimization 
by strengthening the structure to the identified critical response 
and elevating all other design constraints to the level of their 
maximum allowable limitations. Furthermore, all penalties 
values for the best ranked recorded solutions are zeros. That 
demonstrate the competency of the developed algorithm to 
exclude all potential solutions which violates the design 
constraints requirements (Ahmed et al., 2022; Aslam et al., 
2022; Faraj et al., 2022). 
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3.3.1 Optimization of structural sections 

Selected structural sections from the pre-defined variation list 
for different structural elements are assembled to form an individual 
solution for the lateral system of the structure. Variation of the 
sections are shown in Figure 13, Figure 14 and Figure 15 for core- 
wall, columns and beams sections, respectively. It can be noticed that 
during the first rounds of generations, section varied greatly in sizes 
to find an optimum combination that represent the objective 
function. Another observation is that most of the change in 
sections occurs among core wall sections. That is justified by the 
importance of core-wall system in resisting lateral loads in 
comparison to other structural elements (Khan, 2018; Khan 
et al, 2019; Alkhatib et al., 2021b; Khan, 2021; Najm et al., 2021; 
Najm et al., 2022). 


3.3.2 Optimization of structural sections 

Finally, the recorded weight minimization for every generated 
solution is shown in Figure 16 in terms of the concrete volume for 
the lateral structural system (Al-Tayeb et al., 2022; Almeshal et al., 
2022; Emad et al., 2022; Saeed et al., 2022). After generating the first 
population, a huge drop of the structural weight is experienced due 
to the randomness of the first population where optimization 
algorithm easily can search for the optimal solution with 
After that, 
generated through 


minimum sections. refinement of the selected 


solutions is selection procedures and 
evaluation techniques to further optimize the structural until 
convergence is reached and all generated solution beyond that 
represent the same section sizes (Khan, 2017; Najem and 
Ibrahim, 2018). 

Optimal sections weight is then compared to the original design 
to compute the reduction in concrete volume as shown in Table 3. 
The effectiveness of this optimization approach is validated through 
the achievable reduction of lateral structural concrete volume by 
11,204 cubic meter of concrete that reflects a 36.1% of the total 
volume from the original sections. Considering the sustainability 
value in other hand with almost 410 Kg of CO, emissions is 
produced for each cubic meter of concrete, the optimization was 
able to eliminate a total of 4,400 tons of embedded CO, from this 


process. 


4 Conclusion 


This study aimed to extend the research and development of 
structural optimization with focus on irregular tall buildings 
subjected to wind motions. The developed approach consists 
of two major procedures. Firstly, an appropriate computational 
technique to evaluate wind loads on irregular tall building is 
proposed. That is achieved by integrating parametrical model 
into computational fluid dynamic analysis, that is coupled with 
fluid structure interaction algorithm to produce wind forces. The 
generated forces are then used as an input, together with gravity 
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forces, to perform structural optimization. Structural 
optimization is developed in this study by employing genetic 
algorithm with enhanced constraint’s function. Mathematical 
formulation of the enhanced design constraints is developed with 
consideration to lateral displacement, drifts and top acceleration 
for regular and complex motion that may be accounted due to the 
irregularity of building. Finite element analysis is then performed 
for each generated solution to evaluate static and dynamic 
structural responses against the pre-defined objective function. 

Existing case study of 75-story building with irregular twisting 
shape is performed to validate the feasibility of this research work. 
Only steady and incompressible flow applied at (0°) angle of attack 
was implemented in the CFD model to simulate the wind flow on 
the studied building. The developed approaches proved to be 
effective in optimizing the lateral structural system of the studied 
building. That is experimented by reducing the concrete volume by 
35.71% from the original design. Hence, It can be concluded that the 
adopted computational methods contribute significantly to the 
performance of design engineers when dealing with complex and 
irregular structures. 

Future studies shall be extended by employing different CFD 
model including flow type and mesh resolution to improve the 
accuracy of capturing the wind motions on tall buildings. 
Moreover, further examination of different evolutionary 
algorithms, such as Neural Network and Practical Swarm, can 
be done to measure their effectiveness in dealing with structural 
optimization problems. Another suggestion, with the benefit of 
parametrical approach employed in this study, another type of 
optimization can be developed and integrated within the 
proposed framework. That includes aerodynamic optimization 
to mitigate wind motions on the building, and topology 
optimization. To evaluate different lateral structural system for 


wind resistance optimal solution. 
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